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We study the renormalization of a single impurity potential in one-dimensional interacting electron 
systems in the presence of magnetic field. Using the bosonization technique and Bethe ansatz 
solutions, we determine the renormalization group flow diagram for the amplitudes of scattering of 
up- and down-spin electrons by the impurity in a quantum wire at low electron density and in the 
Hubbard model at less than half filling. In the absence of magnetic field the repulsive interactions 
are known to enhance backscattering and make the impurity potential impenetrable in the low- 
energy limit. On the contrary, we show that in a strong magnetic field the interaction may suppress 
the backscattering of majority-spin electrons by the impurity potential in the vicinity of the weak- 
potential fixed point. This implies that in a certain temperature range the impurity becomes almost 
transparent for the majority-spin electrons while it is impenetrable for the minority-spin ones. The 
impurity potential can thus have a strong spin-filtering effect. 

PACS numbers: 71.10.Pm, 73.63.Nm, 71.10.Fd 
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I. INTRODUCTION 

Quantum transport in one-dimensional (ID) electron 
systems has been a subject of great interest for many 
years. In one dimension the interplay between electron- 
electron interaction and residual disorder is a crucial fac- 
tor determining transport properties. It is well known 
that the repulsive interactions between electrons strongly 
enhance the backscattering of electrons by impurities at 
low temperature.^'^ This phenomenon is a manifestation 
of the fact that ID electron systems at low temperatures 
become Tomonaga-Luttinger (TL) liquids, with proper- 
ties very different from those of conventional Fermi liq- 
uids. The renormalizations of the electron scattering by 
impurities have been observed experimentally by measur- 
ing temperature or bias dependence of the conductance 
of quantum wires%^. and carbon nanotubes.-^^^ 

In this paper we explore the TL-liquid renormaliza- 
tions of the potential of a single impurity in the pres- 
ence of a strong magnetic field B. Such a field causes 
significant polarization of electron spins. This polariza- 
tion modifies the low-energy properties of the TL liquid, 
resulting in qualitatively different renormalization group 
(RG) flows of electron backscattering by the impurity. In 
particular, when the energy band is less (more) than half 
filled, repulsive interactions may decrease the backscat- 
tering of majority-spin electrons (holes) in the vicinity of 
the weak- impurity fixed point, while increasing that of 
minority-spin electrons (holes). This can be thought of 
as enhancement of a spin-filtering effect due to interac- 
tions. 

Renormalizations of the electron scattering by an im- 
purity at zero magnetic field has been studied analyti- 
cally in the limiting cases of weak and strong scattererii*^ 
The magnitude of the impurity scattering of spin-A elec- 
trons (A =t, i) is characterized by a small backscattering 



amplitude v\ in the former limit, and by a small ampli- 
tude t\ of tunneling through the impurity potential in 
the latter one. The interactions between electrons give 
rise to power-law renormalizations of these amplitudes at 
low temperatures 
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where T is the temperature and D is the bandwidth. The 
exponents ax and at i? = are determined by the 
electron-electron interactions. The repulsive interactions 
result in = ax < and /3| = /3| > 0. Thus, as the 
temperature T is lowered, the scattering of electrons by 
an impurity is enhanced in both the weak and strong 
impurity limits. 

Quantitative resultsi*^ for the exponents ax and (ix 
have been obtained by using the bosonization technique. 
This method provides a convenient description of the 
low-energy properties of ID electron systems in terms 
of bosonic fields (j)x{x) and 11^(3::) satisfying the commu- 
tation relations 



[cpx{x),nx'{x')]=i5xx'5{x-x') 
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The effective Hamiltonian at £? = has the spin-charge 
separated form H = Hp + with 
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Ha — J Hadx, H. 
Here the fields 
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(Hi ± III )/ \f2 describe excitations of the charge and spin 
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modes, Up and are the velocities of the charge and spin 
excitations, and Kp and K„ are the TL-hquid parame- 
ters. It is known that for repulsive interactions Kp is 
smaller than 1, while Kcr scales to 1 at low energies as 
required by the SU(2) symmetry of the problemiiSi 

The results for the renormalizations of the impu- 
rity potential are obtained by adding to the Hamiltonian 
the perturbations describing the impurity scattering and 
tunneling processes, and then studying their scaling di- 
mensions. At i3 = the exponents are related to the 
TL-liquid parameters as 



Kp + K, 
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Hence, at Kp < 1 and K^^ ~ 1, the exponents a\ is neg- 
ative while f3\ is positive. It is interesting, however, that 
if the condition Ka- = 1 imposed by the SU(2) symmetry 
of the problem is relaxed, the sign of the exponents ax 
and /3a may change depending on the value of Kp < 1. 
For example, for a system with spin anisotropy resulting 
in Ka > I and 2 ~ K^ < Kp < K„/{2Kc, ~ 1), one has 
OA > and /3a > Qt^'^ In this case, the potential of a 
weak impurity becomes an irrelevant perturbation, and 
the strength of the scatterer scales to zero at low tem- 
peratures. On the other hand, a scatterer with strength 
exceeding a certain critical value grows at T ^ 0, as 
indicated by positive exponent /3a- 

This analysis suggests that the renormalizations of im- 
purity potential may lead to a suppression of a weak 
impurity when a magnetic field is applied to break the 
SU(2) symmetry. However, the effect of magnetic field 
cannot be correctly described by using the Hamiltonian 
and Q since the spin-charge separation in the ID 
system is destroyed by the magnetic field liSiii Thus, in 
order to understand the dependence of the exponents a\ 
and (3\ on magnetic field, it is necessary to generalize 
the low-energy Hamiltonian to the case i? ^ 0. Such a 
generalization was accomplished in Ref. 12 where the ef- 
fective Hamiltonian for weakly interacting electrons was 
obtained as a Gaussian model consisting of two indepen- 
dent branches of bosonic excitations. Using the effective 
Hamiltonian, one can investigate the scaling of both the 
weak backscattering due to an impurity, and weak tun- 
neling through the impurity potential. We will see that 
if the magnetic field is sufficiently strong, a spin-filtering 
phenomenon, in which the impurity blocks the transport 
of minority-spin charge carriers and only weakly scatters 
the majority-spin ones, can be realized in a certain regime 
of the RG flow. 

The paper is organized as follows. We discuss the 
bosonization approach to ID interacting electron systems 
in a magnetic field in Sec.m The general form of the ef- 
fective Hamiltonian is presented in Sec. Ill Al In Sec. Ill BI 
we obtain the parameters in the effective Hamiltonian for 
the ID Hubbard model from their Bethe ansatz integral 
equations. In Sec. Ill Gl we derive the effective Hamilto- 
nian for a model describing quantum wires at low elec- 
tron density. In Sec. lIIII we discuss the renormalization of 



impurity potential v\. The scaling dimensions of the im- 
purity potential and tunneling operators are calculated 
in four limiting cases in Sec. IIII Al The RG flow dia- 
grams are discussed in Sec. IIII Bl Section Hvl is devoted 
to summary. 



II. EFFECTIVE HAMILTONIAN 
A. Bosonization approach 

In this section we review the low-energy eflFective the- 
ory of ID interacting electrons in a magnetic field, follow- 
ing and extending the bosonization approach introduced 
by Penc and S61yom>i^ This will be the basis of our anal- 
ysis in the following sections. 

To construct the low-energy theory, we first take the 
noninteracting part of the Hamiltonian and linearize the 
dispersions around the Fermi points -iikp\. The electron- 
field operators are expressed in terms of chiral fields as 



where "i! hx (^'la) is the field of right- (left-) moving elec- 
trons of spin A. Following the standard procedure^*^ the 
chiral fields are bosonized as 



*pa(x) 



where 
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Ka is the Klein factor satisfying {k\,k\i} — 25\x and 
k\ — K^, and a is a short-distance cutoff. The bosonic 
fields obey commutation relations 

[(fiPX (x) , ippx> (y)] = inspSxy sgn{x - y) , 

[fRX (x) , (PLX' {y}] = -iTrSxy ■ (9) 

The Hamiltonian density of noninteracting electrons is 
given in terms of the chiral fields by 



(10) 



where ua > are the velocities of the linearized disper- 
sion of the spin-A branch. 

The interactions between electrons result in two- 
particle scattering processes. In the most general case, 
they can be classified^ into the following four types: 
backward scattering (the gi process), forward scattering 
((72), Umklapp scattering (173), and scattering within one 
branch (54). Among these scattering processes, the gi 
interaction between electrons with opposite spins can be 
discarded since kp^ ^ kpi- The gi interaction between 
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electrons with equal spins is equivalent to the (72 scatter- 
ing. Furthermore, since the electron density is assumed 
to be incommensurate with the lattice, the Umklapp ((73) 
scattering can be ignored. As a result, the most general 
form of the quadratic part of the interaction Hamiltonian 
density reads 



n 



int 



92X- 



+S4A 



dx dx 

d(fpx d(pp\ 
dx dx 
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dx dx 
dippx dippx 



dx dx 



where the coupling constants gix.± are real, and 

T forA-i, 



L for P = R, 
R for P = L. 



(11) 
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Combining Eqs. 1)10(1 and Hll() . we find the total effective 
Hamiltonian density written in the matrix form 



47r 



where ip^ = {(pp-^ , lpli , (fipi, (fLi) 
ii = 



dxip^{x) ii dj:ip{x) 



(13) 
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Mi + .94i 



(14) 



and gix,i_ = gix,±/i2Tr). 

We show in Appendix El that the matrix H can be 
brought to the form 



(15) 



with real vectors ijJpi, satisfying the orthonormal condi- 
tions 



[ijjpi,)^ Cijjptyi — 3p6ppi6yi, 



(16) 



Here the subscript v takes two possible values, which we 
will denote as c and s, parameters Uc and Ug are positive, 
and the matrix C is defined as C* = diag(l, — 1, 1, — 1) 
and accounts for the sign factor sp in the commutation 
relations (O- 

We then introduce chiral fields 



(17) 



satisfying the same commutation relations as the original 
fields tfpx, see Eq. In terms of these new fields the 
Hamiltonian H13|l takes the simple form 

^ = 1^ E M[d.VRAx)Y + [d.^L.{x)f}. (18) 



The positive constants and have the meanings of 
the velocities of the two types of elementary excitations 
of the Hamiltonian Ii. We refer to these excitations as 
the holons and spinous. We then introduce the fields 



1 



{(pLu + PRu), 



H, 



1 



-dxilfiLiy - IfiRu), (19) 



and rewrite the effective Hamiltonian density as 

2" 



27r ^ 

v—c^s 



^2tt2 
TT 11„ 



(20) 



Hence the system of ID interacting electrons in a mag- 
netic field can be described as a two-component TL 
liquid. — We note that at i? = the holon and spinon 
modes reduce to the charge and spin modes in Eqs. 
and respectively. 

The fields (f>i, and H^ are related to the original bosonic 
fields 

<Px^7:(vlx + Vr\), = i^dxiipLX- ifiRx) (21) 

Z ZTT 

used in the Hamiltonian (j3Jl, Q through Eqs. H17|) and 
(|19|l . Due to the parity symmetry of the system, this 
linear relation is simplified to 



He 
H., 



, (22) 



where the real matrix A can be obtained from the vectors 
ujpi,, see Appendix IbI Note that our Hamiltonian H20() 
is determined by 6 parameters: the velocities of holons 
and spinous Uc and Us, as well as the four elements of the 
transformation matrix A. 

An alternative approach to bosonized description of ID 
systems in magnetic field was used in Ref. il4. In their 
theory the Hamiltonian depends on 5 parameters, namely 
the velocities Up and Ua- of the charge and spin modes in 
the absence of the field, the TL-liquid parameters Kp and 
Kcr, and the difference of velocities of spin-t and spin-| 
electrons induced by the field. (In the realistic case of 
spin-independent interactions between electrons, — 1, 
and the number of parameters is further reduced to 4.) 
This simplification of the theory^** occurred because the 
magnetic field dependence of the coupling constants de- 
scribing the electron-electron interactions was neglected. 
We believe the approach of Ref. 14 is therefore inappli- 
cable beyond the regime of weak magnetic field. 



B. One-dimensional Hubbard model 

The effective Hamiltonian (|20|) depends on six param- 
eters, Uc, Us, and the four matrix elements of A. In the 
case of exactly solvable models, these parameters can be 
obtained exactly by solving the integral equations of the 
Bethe ansatziiSiii Here we discuss one such case, namely 
the ID Hubbard model. 
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The original Hamiltonian of the ID Hubbard model 
has the form 



B 



(23) 



where c\ is the creation (annihilation) operator of 

spin-A electron at the Ith site, ni,\ = cj ^^Ci ^, the matrix 
element t accounts for hopping between neighboring lat- 
tice sites, and U is the strength of the on-site repulsion; 
t,U > 0. Throughout this paper, we concentrate on the 
case of less than half-filling, i.e., the total electron den- 
sity n < 1. We will comment on the case of n > 1 in 
Sec.HVI 

The ID Hubbard model allows for exact solution by 
Bcthe ansatz at arbitrary value of the field B. This so- 
lution enables one to obtain not only the velocities Uc 
and Us, but also the asymptotics of various correlation 
functions at large distances. The latter are expressed in 
terms of the so-called dressed charges Z^,yi , which can be 
found exactly by solving integral equations of the Bethe 
ansatz^fliii For example, the dressed charge matrix is 
given by 



Zs 



e 

|e 1/V2 



at i? = 0, and 



Zee Ze. 

Z.ee Zs. 




(24) 



(25) 



in the saturation limit where the electron spins are fully 
polarized. Here the dressed charge ^ is defined in Eq. 
(5.1) in Ref. and takes values in the range 1 < ^ < 
•\/2. As the magnetization increases, the dressed charges 
change continuously between the values in the limiting 
cases. 

By comparing the critical exponents obtained from the 
Bethe ansatz with those from the effective Hamiltonian 
(|20|1 . one can relate the matrix A to the dressed charges 



A 



(An 


A12 ^ 




I A21 


A22 ) 





Ze 
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Zs 
Zr 



Zsc 

-Zs 



(26) 



We will use these results in Sec. IHll 

As a simple model for an impurity potential in the 
Hubbard model, we can take the on-site potential 



VHub = Vb(no,| -I- no,x)- 



(27) 



To analyze low-energy transport properties, we take the 
continuum limit, where the density operator is approx- 
imated as no, A ~ ^'|(0)^';^(0). With the chiral elec- 
tron fields the potential is further reduced to the form 
V = V^ +Vi, where 



Vx = vx 



*L(0)*LA(0) + *L(0)*flA(0)l . (28) 



Here we have kept only the backward scattering terms 
and discarded the forward scattering ones, as the lat- 
ter do not affect the conductance. We have also intro- 
duced spin-dependent backscattering amplitude vx- Fi- 
nally, Eq. H28|l is written in terms of the bosonic fields as 



Vx^-—cos[2(^x{x^Q)]. 
na 



(29) 



We set the amplitudes vx to be positive, which is always 
possible by the transformation (/)a ^ '/>a + const. 



C. Quantum wires at low electron density 

Here we derive the low-energy effective Hamiltonian of 
a quantum wire in the low-density limit, where the ef- 
fective Hamiltonian H2U|) takes a particularly simple form 
similar to Eqs. Q and Q), as we will see below. 

When the electron density in the wire is very low, the 
electron-electron interactions are effectively very strong. 
In the limit of infinitely strong repulsion, electrons can 
never occupy the same position in space and can be 
viewed as distinguishable particles. As a result, the en- 
ergy of the electron system becomes independent of the 
electron spins. At strong but finite interactions, the elec- 
trons in the wire can exchange their positions, and the 
spins of neighboring electrons are weakly coupled to each 
other. The resulting spin dynamics is described by the 
Heisenberg model. 



Ha — J X Si ■ Sl+i. 



(30) 



Hence the Hamiltonian of the wire at zero magnetic field 
B — takes the spin-charge separated form-^^'^^ H = 
Hp + Ha- with the two terms given by Eqs. Q and H3U|I . 

At energy scales below the exchange constant J the 
Hamiltonian (|30|) can be bosonizedf2i& and the form © 
of the Hamiltonian density Tia is recovered. The ad- 
vantage of using the Heisenberg form H3U|I is that the 
magnetic field B can be easily incorporated by adding 
a term —\g\^BBS'^ , where g is Lande factor and /x^ is 
Bohr magneton. The field B polarizes the spins and re- 
sults in finite magnetization.^^ In the following, it will be 
convenient to parametrize the Hamiltonian by a relative 
magnetization m defined as m = (n^ — ni)/{n^ -f n^), 
where n^^i are the densities of electrons with given spin 
components. At m < 1 the Hamiltonian of the Heisen- 
berg model in a magnetic field can be bosonized to the 
form with the velocity Ua and the coupling param- 
eter Ka becoming functions of mi^ As m varies from 
to 1, the velocity Ua{m) changes from ■KJ/[2Tl{n^^ + rij)] 
to zero, and Ka{m) grows from 1 to 2. 

Using the separation of charge and spin variables in 
the form © , (|30|) and above mentioned properties of the 
Heisenberg model, one can conclude that the low-energy 
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Hamiltonian density of strongly interacting electron sys- 
tem in a magnetic field has the form 



n = 



hu„ 



TTKp{Up + niRa 



hua{m) 



'KK„{iin) 



(31) 



The first line of the Hamiltonian H31|) describes the charge 
density excitations (holons) of the electron system. Since 
the coupling of the spins is very weak, the magnetic field 
polarizing the spins does not affect the speed of holons Up. 
The form of the holon part is thus essentially equivalent 
to Eq. 10), with the addition of the term toIIct to the 
momentum density. This correction does not affect the 
dynamics of the holons, as [XIo-, dx(f)p\ — 0. On the other 
hand, the addition of mllo- to the momentum density 
ensures that the holon wave carries the spin current due 
to the finite magnetization m of the ground state. Indeed, 
the equation of motion for the holon wave results in the 



relation (hr, ~ 



between the spin and charge currents. 



The form of the spinon part of the Hamiltonian es- 
sentially reproduces the bosonized Hamiltonian of the 
Heisenberg model at finite magnetization; in particu- 
lar, the dependences Mn-(m) and K„{m) are equivalent 
to those discussed in Ref. The only difference is the 
addition of the term — m dx4'p to the spin density. Due 
to the commutation relation [0p, H^] = 0, this term does 
not change the spin dynamics. However, its presence en- 
sures that in the spinon ground state the spin and charge 
densities are proportional to each other: dx4>a = fTidx4>p- 

The effective Hamiltonian can be easily brought 
to the standard form H20|) . with the matrix A taking the 
form 



TABLE I: Scaling dimensions of the backscattering and tun- 
neling operators V\ and Tx in the four limiting cases I33all - 
(I33dl l. The expressions of xv\ and xt\ are given in Eqs. I40II 
and 1451 . respectively. 





limits 










Ti 


(a) 




< D 


Xv\ 


Xvi 






(b) 


\M\Ati\ 


<^D 






Xti; 


XTl 


(c) 




< D 


1/xt] 








(d) 




< D 




1/xTi 







with the distance between the electronsii^ [Here d is the 
distance from the wire to the nearest gate, and as is 
the effective Bohr radius in the material; as ~ 10 nm in 
GaAs.] In this special case the electrons in the quantum 
wire can be described^ by the Hubbard model in the 
limit of low filling, n — > 0, when the discreteness of the 
lattice can be neglected. In particular, the fact that the 
spin excitations arc those of the Heisenberg model 
corresponds to the well-known property-^ of the Hubbard 
model in the limit U/t ^00. In this limit the param- 
eter Kp takes the value 1/2. We have checked that at 
Kp — 1/2 our result (|32|l for the matrix A coincides with 
the result H26() for the Hubbard model with the dressed 
charges Z^j^i found in the limit C//i — > 00 in Ref. 12^ . 

Similarly to the case of the Hubbard model, the pres- 
ence of an impurity in a quantum wire results in backscat- 
tering of electrons. We will use expressions H28|l and 129|) 
for the perturbations in the Hamiltonian, cf. Refs. (ijy. 



A 



All A12 
A21 A22 



^(1 



K„{m) 



(32) 

In general, the parameter Kp is non-universal. In the 
limit of strong short-range interaction it can be deduced 
from the well known properties of the Hubbard model, 
and one finds Kp = i. For longer range interactions 
one expects Kp < i. On the other hand, the parame- 
ter Kc{m) is the TL-liquid parameter for the Heisenberg 
spin chain in magnetic field, which can be determined ex- 
actly by solving the Bethe ansatz integral equations>i2i2£ 
These results will be used in Sec. IIIII to investigate the 
renormalizations of impurity potential. 

Our discussion in this section assumed arbitrary range 
of interactions between the electrons in the quantum 
wire. In experiments the range of the Coulomb repulsion 
between electrons is usually longer than the distance be- 
tween particles. However, the range of the interactions is 
limited by the presence of metal gates in the vicinity of 
the wire. One can show that at electron densities below 
aBl<P the range of the interactions is short compared 



III. RENORMALIZATION-GROUP ANALYSIS 



In this section we discuss the renormalizations of im- 
purity backscattering amplitude vx and tunneling am- 
plitude t\ using the effective Hamiltonian obtained in 
Sec.m We will consider the following four limiting cases: 



Iwil < D and \vi \ < D, 
|t|| < £1 and |tx| < D, 

< L* and |tj < D, 

< D and \vi \ < D. 



(33a) 
(33b) 
(33c) 
(33d) 



Scaling dimensions of the backscattering and tunneling 
operators in the four cases are summarized in Table 
Evaluating the scaling dimensions quantitatively, we con- 
struct the RG flow diagram in the v^-v^ plane. We will 
see that the RG flow diagram changes drastically when 
a sufficiently strong field is applied. 
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A. Scaling dimensions of backscattering and 
tunneling operators 

1. Weak-potential limit |ux] <C -D 

First, we discuss the renormalization of the impurity 
scattering in the hmit where both V'\ and are weak, 
I^tIj kil ^ To find the scahng dimension of the oper- 
ator H29|l . we calculate the ground-state correlation func- 
tion of V\ using the action for the pure system H2()|l given 

by 



S 



E 



dr I dx 

J-oo 



+ 



1^ (dar:4>uix,' 



,(34) 



where (3 — 1/T and r is the imaginary time. Since V\ 
depends only on </i°(t) = (j)x{x = 0,t), we integrate out 
the fields 4>i,{x,t) a,t x ^ to derive the effective action 
for (cc = 0, r). 



(35) 



Here we have introduced the Fourier transform 
1 rK 



0^(a; = O,r)dT (36) 



with LUn = 27rri//9. Using the effective action, the 
imaginary-time correlation function is calculated as 



/e2«<^'A(^)e-2»0A(o)\ 



I?0j!l?.^2e-^°+2^[^"(^)-*"(")l, (37) 



where 



(38) 



is the partition function. A straightforward calculation 
of Eq. (PTjl with Eqs. ll^ and gives the correlation 
functions in the limit /3 ^ oo, 

(V^a(t)Fa(O)) (39) 

with the scaling dimensions xyx given by 

xy| = Aj^+Ali, (40a) 

(40b) 



XVl — + ^422- 



The exponents a\ in Eq. Q in this limit are given by 
= ^vx - 1- (41) 



2. Weak-tunneling limit \tf\, \ti \ <^ D 

The renormalizations of tunneling through a strong 
impurity potential were studied using several different 
approachesii*2i2^ The discussion given below uses the 
method of Ref. IH 

Let us consider the tunneling of a spin-]" electron 
through the impurity potential at a; = 0. Since the poten- 
tial amplitudes v\ (A =1,1) are assumed to be very large, 
the fields (p*^ are pinned at the minima of the potential V\ 
[Eq. i.e., 0° = nix, where l\ are integers. The tun- 

neling of a spin-t electron through the potential barrier 
is equivalent to a sudden change of (j)'^ between neighbor- 
ing minima, say, from (f)'^ — to (j)'^ = tt, i.e., a jump in 
0| by TT. Let us denote the operator for this tmmcling 

process by T|. The correlation function {T^{t)T^{Q)) is 
then obtained from 



(T/(r)T^(0))cxexp(-5o) 



(42) 



where we have ignored small fiuctuations of 0^ around 
the potential minima. [Here 9(t) is the unit step hmc- 
tion.] From Eq. H22(l we substitute 



(43) 



into exp(— 5o) to find 

(rt(r)r^(O)) «r-2--T, (44) 
where the scaling dimension is 

XT, = [{A-'W + [{A-')uf = ^f^- (45a) 



(detA)2 



Similarly we have 



xTi^[{A-^ur+[{A-^h2r^^^^. 

(det A)^ 



(45b) 



The exponents for the tunneling amplitude t\ in Eq. 
in this limit are 



/3f =XTA-1. 



(46) 



3. Asymmetric limits \vi\,\ti \ -C D and \t^\,\vi \ <^ D 

Next we consider the asymmetric limit where the po- 
tential scattering is weak for the spin-| electrons but 
strong for the spin-| electrons. Although such an ex- 
tremely spin-sclcctive scattering is not likely to be re- 
alized with a bare impurity potential, we will see in 
Sec. lHI Bl that this is indeed realized for some RG trajec- 
tory if a magnetic field applied is sufficiently large. 

The scaling dimension of the potential in the limit 
ktl' I'^il ^ -D can be found in a similar way to the weak 
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potential limit discussed in Sec. IIII A II The only differ- 
ence is that in the present case (piix — 0) is pinned at a 
potential minimum ttI^ by the strong impurity potential 
Vi . The asymptotic form of the ground-state correlation 
function is then obtained as 



i>°=o 



(47) 



with XTf given by Eq. ()45a|l . Therefore the scaling di- 
mension of V-\ in this limit is I/xt] -, and the exponent 
a-f is given by 



1 



1. 



(48) 



Similarly, the scaling dimension of Vi in the limit 
KtI' l"*^!! ^ ^ is found to be 1/xti, resulting in the ex- 
ponent a'"^^ — 1/xti ~ 1- 

The scaling of the tunneling operator T| in the limit 
KtI' l^il ^ ^ studied in a similar manner to that 

in Sec. Ill A2. In the present case, however, the potential 
Vi is weak, and the field 0° can fluctuate almost freely. 
Therefore, to find the scaling dimension in lowest order in 
vi, we first integrate out in 5*0 to obtain the effective 
action for 0°, into which wc substitute 0|(t') = TTd{T' — 
t). This yields 



(Tt(r)r^(O)) cx j V4>le 

-2/a;vT 



!,0(T')=7re(T'-r) 



CX r 



(49) 



with xv] given by Eq. (|40a|l . We thus conclude that 
the scaling dimension of T| in the limit Ij33d|) is l/xyf, 

and the exponent (3^^ = l/xv\ 
ing dimension of Tj in the limit 
exponent /jj*^' — l/xv\, ~ 1- 



1 



Similarly, the scal- 
c|) is Ijxvi, and the 



B. RG flow diagram 
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FIG. 1: Magnetization dependence of scaling dimensions xv\ 
and xt\, computed from the Bethe ansatz integral equations, 
(a) for the Hubbard chain with U/t = 4 and n = 0.5 and (b) 
for the quantum wire with Kp = 1/4 in the low-density limit. 



In the preceding sections we have found that the scal- 
ing dimensions of the backscattering and tunneling op- 
erators in the four limits H33a() - (|33d|l are given in terms 
of the matrix elements of A. We thus need to compute 
these matrix elements to determine whether the perturb- 
ing operators in each limit are relevant or irrelevant. As 
we have seen in Sec. lUl this can be achieved for electron 
systems with short-range interactions, i.e., the Hubbard 
chain at less than half-filling, and the quantum wire in 
the low electron-density limit. The dressed charges Z^i,/ 
for the former case and the parameter (m) for the lat- 
ter one can be calculated as functions of m by solving the 
corresponding Bethe ansatz integral equations. We have 
solved the integral equations numerically, and the results 



are discussed below. The readers who are interested in 
the details of the Bethe ansatz analysis should refer to 
Refs. ^3 and ^ for the Hubbard chain and Ref. |23 for 
the Heisenberg chain. 

In Fig. n we plot the scaling dimensions xyx and xtx 
obtained numerically for the Hubbard chain at n < 1 and 
the low-density quantum wire as functions of the relative 
magnetization m. The two systems exhibit qualitatively 
the same behavior. At zero magnetization m = 0, we see 
that xv^ — xvi < 1 and xt]- — xti > 1. This is in ac- 
cordance with the well-known scaling behavior in SU(2) 

symmetric systems, aj'*'' = a^'*'' < and — p'"^^ > 0. 
The dimensions xyf and xyi increase with m. In the 
limit of full spin polarization m — > 1 the dimension 
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FIG. 2: Schematic RG flow diagram for (a) m — 0, (b) < m < nic, and (c) m > rric. 



xv\ reaches a certain value y greater than 1, whereas 
XVI ~^ 1- (I"^ Hubbard model xti and xt\ also ap- 
proach y and 1 as m — > 1.) The most important point 
here is that the dimension xv\ exceeds 1 for m larger than 
certain critical magnetization 771^. This implies that the 

exponents a^-^^"^ and /jj'^'' change their sign at m = nic, 
and the direction of the RG flows reverses. The signifi- 
cance of this effect can be quantified by the value of y, 
which is given byy=l + {l — -I arctan[4t sin(7rn)/{/]} 
for the Hubbard chain at n < 1 and y = 1 + 2Kp for the 
quantum wire in the low electron-density limit. On the 
other hand, the dimensions xt] and xti are larger than 
1 at any m, indicating that the scaling of the operators 
related to them does not change qualitatively between 
m > rric and m < nic- We have checked that the de- 
pendence of the scaling dimensions on the magnetization 
m remains qualitatively the same regardless of the inter- 
action strength U/t and the electron density n < 1 in 
the Hubbard model, or the exact value of Kp < 1/2 for 
low-density quantum wires. 

From the magnetization dependence of the scaling di- 
mensions xv\ and xt\ discussed above, we can deduce 
the RG flow diagram of {v^ ,vi) as shown in Fig.|21 When 
the magnetization is small, m < mc, the backscattering 
of electrons by an impurity is enhanced by repulsive inter- 
actions. As a result, the RG trajectories go directly to the 
strong-backscattering fixed point = (cx),cx)) [Fig. 

|2Ia) and (b)]. The magnetic field bends the RG trajec- 
tories upward but does not change the essential features 
of the renormalization flow. 

On the other hand, when the applied field is suffi- 
ciently strong to achieve m > rric, the backscattering 
operator V| of spin-l electrons becomes irrelevant in the 
weak-potential limit. That is, the electron-electron in- 
teractions suppress the backscattering of majority-spin 
electrons by a weak impurity. As a result, the RG trajec- 
tories in the vicinity of = (0,0) flow toward the 
line = [Fig. [2Ic)]. Thus, if the bare backscattering 
amplitudes are not too large, they are renormalized to- 



ward the asymmetric limit (uj, v^) = (0, 00). This means 
that, in a certain regime of RG transformations or equiv- 
alently, at certain energy and temperature range, a sit- 
uation is realized where the impurity potential becomes 
almost transparent for majority-spin electrons but almost 
impenetrable for minority-spin electrons. However, since 
the asymmetric fixed point {vi,vi) = (0, cxd) is unstable, 
the potentials are eventually renormalized to the strong- 
scattering fixed point (w^, vi) = (cx), 00), with decreasing 
energy scale or temperature T — > 0. In this sense, a 
weak impurity potential can have a spin-filtering effect 
generating a spin-polarized current. 



IV. SUMMARY 

In this paper we have studied the effect of magnetic 
field on the RG flow of a single impurity potential in 
ID interacting electron systems. Within the Abelian 
bosonization theory, low-energy physics of the system 
is described as a Gaussian model with two independent 
modes of bosonic excitations. The coupling parameters 
Aij are obtained from the Bethe ansatz for the ID Hub- 
bard model at less than half-filling as well as for a quan- 
tum wire at the low electron density limit. Using these re- 
sults, we have evaluated the scaling dimensions of the im- 
purity potential and tunneling operators, and determined 
the RG flow of the potential amplitudes near the fixed 
points. We have found that the magnetic field can cause a 
drastic change in the RG flow diagram. While in a weak 
field the repulsive interactions always enhance electron 
backscattering by impurities, in a sufficiently strong field 
this effect is reduced, and the backscattering of majority- 
spin electrons by a weak-impurity potential may even be 
suppressed by the interactions. This means that if the 
amplitude of the bare potential is small, a spin-filtering 
phenomenon, in which only the majority-spin electrons 
can transmit through the renormalized potential, can be 
realized in a certain temperature regime in the RG flow. 
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So far we have assumed that the band is less than half 
filled, n < 1. The results obtained at n < 1 can be easily 
translated to the case n > 1 by performing the particle- 
hole transformation, nx —^ 1 — nx. Suppose that the 
band is more than half-filled and partially spin polarized, 
< m < 1. Then, it is the majority-spin holes, i.e., 
the spin- 1 holes, whose backscattering is suppressed in 
strong magnetic field. In other words, transport of charge 
carriers of majority-spin is enhanced by interactions. 

The renormalizations of the impurity potential in the 
presence of magnetic field have also been considered in 
a recent preprint i^iiS^ Although the authors of Ref. |^ 
also found the regime in which the weak backscatter- 
ing of the majority-spin electrons is suppressed by the 
interactions, their results differ significantly from ours. 
Most importantly, in Ref. 2^ this interesting regime oc- 
curs either in the presence of spin-dependent interactions 
between electrons, or when the interactions are attrac- 
tive. Both of these regimes are unlikely to be realized 
in realistic experiments. In contrast to our work, the 
theorji24 was based upon the treatment of ID electron 
systems in magnetic field developed in Ref. ^3- As we 
mentioned in Sec. Ill Al the latter approach assumes weak 
magnetic field, i.e., small relative magnetization m <C 1. 
In contrast, our effect of suppression of weak impurities 
by interactions occurs at sufficiently strong field, when 
iTi > iTic ^ 0.2 (see Fig. ^ and is expected in the realis- 
tic case of spin-independent repulsive interactions. 
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APPENDIX A: GENERALIZED EIGENVALUE 
PROBLEM OF HC 

In this Appendix we show how the matrix (|14|) can be 
brought to the diagonal form (|I5|I . Suppose that ujj and 
T]j are the right and left eigenvectors of HC, respectively, 
i.e., 

HCujj — UjLOj, (Al) 
{ri.fnC = u,{ri,f. (A2) 

Then, the matrix Ti.C is given by 

7i:(7 = ^a;,u,(r,,)^, (A3) 
j 



where u)j and rjj obey the biorthogonal condition, 

iVjf^r^S,,,. (A4) 

Since both H and C are symmetric, i.e., (Ti)"^ ~ Ti and 
{CY = C*, Eq. (|A2(I can be rewritten as 

Hence, rjj satisfies the relation 

ncC'nj = c^nc^-nj = ccH'nj = u^c-nj, (A5) 

where we used = I. This means that Crjj is a right 
eigenvector of "HC and proportional to iVj , 

^3=CjCr]j, (A6) 

where Cj is a constant. Furthermore, assuming that the 
Hamiltonian matrix H is positive definiter& it follows 
that Cj and Uj have the same sign, 

= UjCj > 0. (A7) 

We can therefore relate Tjj to Uj as 

Uj = sgn{uj)CT]j (A8) 

without any loss of generality. Using Eqs. ljA3|l and l|A8p . 
one finds that the right eigenvectors u)j satisfy Eq. H15|l . 

3 

The vectors ujj satisfy certain orthonormal conditions, 
which can be derived as follows. Suppose that P is a ma- 
trix of parity transformation exchanging the right- and 
left-moving fields. Due to the parity symmetry of the 
Hamiltonian, the vectors and Puij have the follow- 
ing property: if ijjj is a right eigenvector of TiC with an 
eigenvalue Uj, i.e., T-LCujj — UjiVj, then Pu}j is another 
right eigenvector, with eigenvalue — Uj, 

nCPuij = -PnCijJj = -UjPijJj, (A9) 

where we used the relations PHP — H and PCP = —C. 
Hence, we can classify the vectors u}j into two pairs, ujpc 
and u}ps, where the right and left movers in each pair 
are related as lj^^, — PijJpy. Using these results and 
Eqs. ljA4p and (jA8|) . we find that the vectors ujpi, obey 
the orthonormal conditions Hlt)|l . 
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APPENDIX B: MATRIX A 



where the matrix A± is given by 



Here we outline the diagonahzation procedure trans- 
forming the Hamiltonian H13|l to the form H20|l and ex- 
press matrix A in terms of the velocities and the 
coupling constants giA,_L- 

We denote the elements of the vectors ojpi, as 



V Oui j 



(Bl) 



Due to the parity symmetry of these vectors, the fields 
(^A and Ha do not mix with each other under the trans- 
formation. The relation between (^^ (n,y) and (^x (IIa) 
has the form, 




= A- 




(B2) 



A, 



(B3) 



with A±^A = o-v\^b^\. Using the orthonormal conditions 
Eq. (|16|l . one can find that the matrices A± are related 
to each other as {A^)^ = (A^)^^. We thereby arrive at 
Eq. by identifying A = A^. 

From the eigenvalue problem Eq. HA1|I , the linear equa- 
tion system for A±^x is obtained as 



/ 



V 54-L - 52_L + 941 - 921 



"T + 54T + 52T 54_L + .g2_L 

54_L + 92± + 941 + 921 







\ 


( ^+.T ^ 




( A+u^ \ 
















A-^X 


/ 


\ A-H ) 







(B4) 



r 



By solving this eigenvalue problem, we obtain analytical 
expressions for the renormalized velocities u,y^ 



and the elements of A 



±1 



A-,i 



K - Pi 



A 



q 



2 \ -1 1 



K-Pi 



uZ-P] (^^ul_-p£\ ^ 1 



K - Pi 



A 



(B5) 
(B6) 



(M| + 94] + 92]){ul -Pl) + {g4± + 92±)Q 



(B7) 
(B8) 

(B9) 

(BIO) 



where 



9l± 



92±, 



P\ = (WA + .94a)^ - .92 A 

9 = (""T + .94T + .92t)(34± - .92±) 

+ ("1 + .941 - .92i)(54± + g2±), 

r = [u-^ + 94] - .92t)(54± + .92±) 

+ + 941 + 92l){94± - 92±)- 



(Bll) 
(B12) 
(B13) 



Substituting Eqs. (6.12)-(6.17) of Ref . 12 into Eqs. (|B5|) - 
(|B13|) . we obtain the expression of matrix A for the Hub- 
bard model [Eq. it^ ]. 
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